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Abstract 

We have compiled some sectors of works by L. Zamick and collabora- 
tors which we hope will be of interest to the reader. Some of the results 
which we, at least, have found amusing often popped up unexpectedly 
in what were apparent routine calculations. We feel that the sections we 
have made can be of importance in farthing the field but our main criteria 
are that the results will have some educational value and that they add a 
bit of zest to our field. 
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Abstract 

A simple formula is obtained, which relates the nuclear compressibility 
to the binding energy per particle, the mean kinetic energy per particle, 
and the power of the density appearing in the two body interaction. 
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In this work we wish to calculate the compressibility of a nucleus, assuming an 
interaction which compressibility of a nucleus, assuming an interaction which 
correctly reproduces the total binding energy and the mean square radius. Such 
that an interaction requires at least two parameters in order to obtain these two 
properties, but we shall use a three-parameter interaction 

-aS (r) + 1P a (R) 6 (r) 

This a familiar density dependent interaction, not unsimilar to the one chosen 
by Moszkowski [1], except that he fixes the power of the density, a, to a value 
of 2/3 and he includes finite range corrections. Vautherim and Brink [2] have 
cr=l. 

Here, on the contrary, the parameter a is kept as a variable, for we wish to 
study what influence it has on the nuclear compressibility. 

We assume that we can restrict ourselves to harmonic-oscillator wave func- 
tions, characterized by only one parameter, the oscillator length b. In that case 
we can by dimensional analysis, write the expression for the total binding energy 
as 

E = -E B = A/b 2 + B/b 3 + C*/6 (3+3cr) 

where A, B, and C are independent of b. The first term is the kinetic energy 
and is given by the oscilator formula 

A/b 2 = (h/ (2mb 2 )) J2 2n + l + 3 / 2 
Since we are at equilibrium 



dE n 



= -2A/b 3 - W/b 4 - (3 + 3cr) C/b^ +3 ^ 

We can regard the above as two equations in the two unknowns B and C. We 
obtain 

B = -b 3 (1 + a) E B jo - (T) b 3 (I + 3cr) /3a 

C = b^+^EB/a + (T) 6( 3 + 3CT ) /3c7 

evaluated at the equilibrium value of b and Eb ■ 
The nuclear compressibility K, is defined as 

We thus obtain an expression for K totally in terms of the binding energy per 
particle and mean kinetic energy of a particle (both positive quantities) 

K = {1/sA) [(T) + qE B + a (3 (T) + qE B )} 

This is our main result. Note that the compressibility is linear in the power 
of the density a and increases with increasing a. 



2 



The above is an obvious explanation of Bethe's observation [3] that the Vau- 
therim Brink interaction [2] (cr=l), yields a considerably larger compressibility 
then the interaction he used. 
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Abstract 

A simple identity is used to show that when harmonic oscillator wave 
functions are used, Inghs formula for the mass parameter of the breathing 
mode state yields the classical result. 

We note the following relationship obey by the operator r 2 (square radius) when 
Slater determinants with harmonic oscillator radial wave functions are used 



El(°ll> 2 h> ^ 2 (o|I> 2 l 



In the above b is the oscillator length parameter, |0) is the ground state of a 
closed shell and |n) is an excited state. We need only consider excited states of 



the type 



.7=0 



that is one-particle-one hole states which differ by a node 



and couple to J=0. Working out both sides we get 

(si + 1) (n + 1) (n + 1 + 3/2) = ^ (21 + 1) (2n + I + 3/2) 



n = 0,1,2... 

On the left we sum only over those occupied states nlj which are such that 
the state n+1, lj is above the Fermi sea; on the right we sum over all occupied 
states. 

The above can be derived by noting that the harmonic oscillator spectrum is 
such that (E n — E D ), the energy difference, is constant —2fiuj. Hence the energy 
independent sum rule for E(0) transitions is proportional to the energy weighted 
sum rule. 

As noted by Werntz and Uberalll [2], if one considers classically a sphere 
of mass M and root mean square radius r\ which undergoes breathing mode 
oscillations then the kinetic energy is 1/2 Mi] 2 . We wish to show that the 
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same result is obtained in the Inglis model [1] using harmonic oscillator wave 
functions. 

Let us first define an auxiliary mass parameter 



M b = 2ft 2 



d_ 



I (Ef Ei 



A harmonic oscillator wave function has the structure N(b)f(r/b), where N is 
the normalization, and b is the oscillator length. We have 

d dN df 

The first term does not contribute because of the orthogonality of |n) and |0). 
Hence we get 

d 



M b =(2ft 2 ) /6 2 ]T 



' dri 







The quantity can be written as a commutator involving the shell model 
Hamiltonian 

T ^t = (- m / 2h2 ) [H sm ,r 2 ]-3/2 

he last term does not contribute. We also can make the simplification (n [H sm , r 2 ] 0) = 
(E f - E t ) (n \r 2 \0). Finally but replacing (Ef - E t ) by 2haj = (2ft 2 ) / (mb 2 ) 
we get 



M b = 2ft 2 



d_ 
db 



/b 4 



But we are interested in the mass parameter M rms , in which the root mean 
square radius r\ is the dynamical variable. We get 



M r 



Hence, we get 



' Mb=(b 2 Mb)/((v\£r 2 \v)/A) 

'E 6 |(«|E- 2 |o)| 2 l 



M rmA = mA ■ 



b 2 (o|E^ 2 |o) 



But from the relationship at the beginning of this paper the quantity in brackets 
is one. Hence M rms = mA. 

For completeness we note the well known equation of motion for the breath- 
ing mode state 

l -M Trns if + X -AK ( V - Vo ) 2 = E 

where K is the nuclear compressibility. The energy of the breathing mode state 
is 

(Mb - [AKI (M rmsV 2 )) 1/2 = y/K/ (mr, 2 ) 
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One often takes rfc = ^r 2 A 2 / 3 with r Q =1.2 fm. I would like to thank J. 
Speth and A. Fassler for their interest and hospitality. 
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Abstract 

The isospin mixing of the lowest two 1* states in 12C is here consid- 
ered. A nuclear charge symmetry breaking interaction (CSB) which was 
previously fitted to help explain the Nolen-Schiffer anomaly, goes in the 
right direction for this problem, but is not as important as in the N-S prob- 
lem. An important effect is mixing due to the one body Coulomb field. 
This is shown to be approximately proportional to (r 2 ) p ^ — (r 2 ) p , 

a quantity which is very sensitive to the details of the one body nuclear 
potential. 

The total off-diagonal matrix elements are 170 keV with the DME and 92 keV 
with the SKII interaction. The DME can reproduce the experimental value 
obtained by Lind et al., while the SKII interaction reproduces the one obtained 
by Adelberger et al. 

Why are the DME and SKII interactions so different? Looking in table 2 
we note that they differ by 62 keV in the off-diagonal matrix element of one 
body coulomb field. Thus the quantity which is conceptually the simplest thing 
entering into this calculation is giving the largest ambiguity. 

By following an approximate procedure of Lane Martorelll and Zamick [10], 
we can relate the matrix element of the Coulomb potential to the square radius 
of a shell model orbit. We approximate Uc by its value in the interior of a 
uniform sphere 



Thus 

(T=l\U c \T = Q) = ^r\(r 2 ) -(r 2 ) 1 (7) 

We now see why we are getting large differences between DME and SKII in- 
teraction. Whereas for thep 3 / 2 proton orbit the difference(r 2 ) £)MS — ( r2 ) SKII — 
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0.126/to 2 or 1.7%, the corresponding difference for the P1/2 proton orbit is much 
bigger, 0.962 fm 2 or 10%. 

Why is there such a large difference between DME and SKII? A possible 
explanation is that the pi/ 2 orbit is bound by only 5.74 MeV for DME, but for 
SKII interaction the value is 7.15 MeV. This certainly goes in the right direction. 

We examine this possibility in more detail by making the following list. 



DefineA = 



( r 2\ _ ( r 2\ 

X 1 Pl/2 X 'PS/2 



proton 



A[/m 2 ] (T = l\U c \T = 0) = ^[keV] 

DME 2.172 178 
WS\DME] 1.844 151 
SKII 1.336 109 
WS[SKII] 1.347 109 

In the above WS[DME] stands for the Woods-Saxon parameters which give 
the same single particle energies as DME etc. Note that A is the same for the 
SKII interaction as for WS[SKII]. However, there is a discrepancy between DME 
WS [DME] such that when all is said and done the DME will give a larger off- 
diagonal matrix element. Perhaps the most interesting comment from the above 
list is that we get a change of 45 keV in going from WS [DME] to WS [SKII]. 
This is, by simply demanding the P1/2 orbit be bound by 5.74 MeV rather than 
7.15 MeV, we are getting a change which is starting to be comparable to the 
difference in the values which have been cited by the two experimental groups. 

In conclusion we note that our results are not inconsistent with the currently 
available experimental data; however, those data have large error bars and there 
is some disagreement between two groups. The phenomenological charge sym- 
metry breaking interaction which helps explain the Nolen-Schiffer anomaly also 
helps to enhance the isospin mixing matrix element. The results depend sensi- 
tively upon the square radii of the single particle orbits. Indeed two different 
Hartree-Fock models give different results for these quantities, and it is hard at 
present to favor one model over the other. Further work along these lines i.e. 
learning what are the correct single particle energies for the shell model orbits, 
is suggested. 
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Abstract 

Isoscalar monopole, quadrupole, and octopole states are calculated in 
closed shell nuclei. A delta interaction is used, its strength determined by 
the multipole condition, namely that the mean single-particle-single-holc 
potential energy difference is equal to the corresponding kinetic energy 
difference. This interaction is used to obtain particle-hole matrix elements 
appropriate for a random-phase approximation calculation. The strengths 
given by different multipole conditions are different but (except L=0) they 
appear to approach each other as the mass number A becomes large. 
The monopole mode has already collapsed somewhat before the strength 
implied by the quadrupole condition is reached. Some of the calculations 
were repeated using a zero-range Skyrme interaction. We observe a very 
high degeneracy in our calculation which we are able to explain in terms 
of L-S coupling and the fact that the particle-hole matrix element of a 
delta interaction is of the form J fp A H A (r) fp B H B (r) r 2 dr. 

IV. EXPLANATION OF HIGH DEGENERACY IN L-S 

Coupling Although a j-j coupling basis was used for these calculations, it should 
be noted that the interaction we use is central, and we are not introducing a 
one body spin orbit interaction. This means that the quantum numbers L and 
S are good. This manifests itself in the fact that most of the states that we 
obtain in an RPA diagonalization have vanishing B(EL)'s to the ground state. 
These correspond to states with S^O. 

Furthermore, if we choose the single particle energies to be those of a har- 
monic oscillator (this was described in more detail earlier as SPII), we find 
that even fewer states than are permitted by L-S coupling have finite B(EL) 
transition to ground states. 



7 



Furthermore, we observe that these states which have vanishing B(EL)'s have 
eigenvalues which are some multiples of hw. For the quadrupole states there is a 
high degeneracy of 2fouj; in the octopole case at lftu and 3fku. Take for example 
L=2 + states in w O. We find that there are 4 states which are degenerate at 
exactly 2hw. These states have zero B(E2) strength to the ground state. This 
degeneracy can be partly explained by looking at the expression for the particle- 
hole interaction in L-S coupling which is given in Appendix A. The particle-hole 
states forming the 2 + states in le O can be divided into 3 classes. 
Class{\) L^2 S = l (0p- 1 0/)i = 3 (Op^lp) L = 1 
Class(2) L = 2 S = (Os^Od) (0p~ 1 0f) (Op^lp) 
Class(3) L = 2 S = l (Os^Od) (Op^O/) (Op^lp) 
The 3j symbol 

lh Ip L 


vanishes unless lh + l p + L is even. Hence all the matrix elements in class 1 
will vanish. This explains two of the four states at 2huj. 

We next note that for every matrix element in class 2 that vanishes there 
must be a corresponding one in class3 that vanishes. This is because the en- 
tire spin dependence is contained in the factor(l — 45s, St,o +^ x ^s,o —2xSt,o )• 
Since we are setting x=0, this factor will be 1 for S=l and -3 for S=0. The 
particle-hole matrix elements in class 2 are -3 times those in class 3. 

Thus two of the four degenerate states are from class 1, one from class 2, 
and one from class 3. The ones from class 2 and class 3 are expected to have the 
same radial structure, differing only in the spins. Only class 2 states can have 
finite B(E2)'s. This is because the E2 operator ^2j 2 y2 has no spin dependence 
and therefore cannot connect spin one to spin zero. Thus, in this example we are 
left with one nontrivial degeneracy (class 2) to explain, as well as the vanishing 
B(E2). 

To explain this, we note that the L-S coupling expression for the particle-hole 
matrix element (given in Appendix A) has the structure 

( [A'" 1 ] LST V [ph- 1 ] LST ) = J f' k (r) h (r) r*dr 

where from here on we use the symbol k to designate ph. 

Let us denote the class 2 state at 2fkJ by -0 and expand it in terms of 
particle-hole components (we limit ourselves to TDA states in this discussion). 

^ = ^2 z k i fc ) 

We Expect 

This can be achieved by demanding 

X] h ( r ) Zk = 

k 
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Let us first consider the schematic approximation in which the integral j R ni i 1 (r) R n2 h ( r ) Rush ( r ) RmU ( r ) 
is replaced by a constant. 4 In this case the particle hole matrix element has the 
structure(k'Vphk) = g' k gk, where g goes not depend on r. We can regard g k as el- 
ements of a vector of dimension D, where D is the number of particle-hole states 
for a given LST configuration, e.g., in the above example, 16 L=2, S=0, T=0 
we have D=3. The quantities z k are also elements of a vector of dimension D, 
and the condition^ fc z k gk — simply means that the vector {z} is orthogonal 
to the vector {g}. 

There are clearly (D-l) vectors {z} which are orthogonal to {</}. Each of 
these (D-l) vectors will have a vanishing particle-hole matrix element. Thus 
the class 2 degeneracy at 2huj will be (D-l). This will also be for the class 3 
degeneracy. Thus for L=2, S-0, and T=0 states in 16 we expect a twofold 
degeneracy from class 2. This was confirmed by a calculation. We obtain the 
same degeneracy in class 3. 

We now consider the case where /& (r) is not approximated by a constant. 
We note that (r) is proportional to the product of two harmonic oscillator 
radial wave functions R n (r), Rh (r). Using the Variable x=r/b we note that 
the product can be written in the form of an exponential times a polynomial 

N max (ph) 

f ph {r) = e- x2 X° Y, a n(ph)X 2n 

N mzn (ph) 

where er=0 or 1. Let iV mi „=mimmum of all N min {ph) and let N max = maxi- 
mum of all N max (ph). To ensure that J2kf k ( r ) z t - " we demand that each 
coefficient oLY 2n+cr vanishes. That is, °« (&) z k = 0. This leads to A con- 
ditions where A = N max — N m i n + 1. We now have A vectors, the nth one 
of which is[a„ (1) , a n (2) , ...a n (D)]. The vector {z (k)} has to be orthogonal 
to all of these. There are clearly (D — A) such vectors {z (k)}. Hence class 2 
degeneracy is (D — A). 

In our example (L=2, S=0, and T=0 in 16 0) the polynomial has terms 
in x 2 and x 4 . Thus N min =l, N max =2, and hence A=2. Thus the degeneracy 
(D — A) is one. 

In a major shell n, 1, N (N=2n+ll) the lowest power that appears in a 
polynomial is x 1 ; the highest power is x N . Thus if the particle is in the shell 
n p , l p , N p , and the hole is the shell rih, lh, Nh, then we have 

N max {ph) = ^(N P + N h ) if N p + N h is even 

N max {ph) = ^(N p + N h - 1) if N p + N h is odd 
N min (ph) = i(Z p + l h ) if l p + k is even 
N m in{ph) = \{lp + lh~ 1) if l P + k is odd 
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Just to give another example, consider the L=2 and T=0 states in °Ca. 
There are seventeen states in all, five in class 1, six in class 2, and six in class 3. 

The entire degeneracy at 2fko for the schematic model 15. Five of the states 
are from class 1. Since D is 6 for class 2 the degeneracy here is D-l=5. It is 
also 5 for class 3. 

For the delta interaction the entire degeneracy at 2hw is eleven. Five still 
come from class 1, leaving 3 from class 2 and 3 from class 3. We find that the 
quantity A is equal to three. The class 2 degeneracy is therefore D — A = 3, as 
expected. 

It is very easy to show that the B(EL) to ground is zero for these degenerate 
states. 

The B(EL) is proportional to J r L J2k fk ( r ) Zkdr. Since the integrand is 
zero, the integral will also be zero. The above is not true when we approximate 
the radial integral by a constant. 

It should be emphasized that while part of the above argument involved spin 
isospin symmetry (as discussed many years ago by de Shalit and Walecka 5 ), the 
crucial part involving class 2 degeneracies did not. 

It is also worthwhile noting that for L-S coupling the expression for the 
particle-hole matrix element the spin-isospin factor is very simple(l — 4(55, o +2x5s l0 —2x8t, )■ 
The values of this factor for the four different modes are -3 for S=0, T=0; (T2x) 
for S=l, T=0; (l+2x) for S=0, T=l; and 1 for S=l, T=l. 
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Abstract 

It is noted that in a single j shell / 7 / 2 calculations in the titanium 
isotopes one obtains Ml rates of about one single particle unit in strength 
to states at about 4 MeV excitation. This is systematically consistent 
with the recently discovered collective excitations in 156 Gd, and neigh- 
boring nuclei. In the single j shell case, through, the spin and orbital 
contributions to the Ml matrix element are nearly the same. 
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Recently a low lying magnetic dipole excitation has been discovered in 156 Gd at 
an excitation energy of 3.075 MeV with a strengths (Ml) t - (1.3 ± 0.2) 
The experiment involving high resolution inelastic electron scattering with the 
Darmstadt linear accelerator. This was followed by the discovery, by the Darm- 
stadt group, of low lying 1+ states in other deformed nuclei such as 154 Sm, 
154 Gd, 164 Dy, 168 Er, and 174 Yb. 

The purpose of this Brief Report is to point 1 out that the presence of An Ml 
state with an excitation energy of 3-4.5 MeV and the order of 1 single particle 
unit of strength may be a more widespread phenomenon. We demonstrate 
this by calculating Ml rates to the lowest 1 + states of the even-even titanium 
isotopes, using the wave function of McCullen, Bayman, and Zamick (MBZ) 14, 15 
and Ginocchio and French. 16 

In a single j shell one can replace the magnetic dipole operator /i by (g p L p + g n L n 
where L p is the angular momentum operator for the protons and L n for the neu- 
trons. We use the same quenched g p and g n as were used originally, 34 g p — 1.50 
and g n — 0.39. The Ml transition must be proportional to (g p — g n ) 2 , for oth- 
erwise, the total angular momentum operator J — L v — L n would be able to 
induce an Ml transition for J-0 to J=l. This is clearly impossible. 

The MBZ wave functions 14 ' 15 for the titanium isotopes are of the form 



* J = J2 D J (L p ,L n ,v)Ufi /2 ) L J (# /2 )* 



L~,L n ,v 



where D J (L p , L n ) is the probability amplitude that two protons couple to L p 
and n neutrons couple to L n . The parameter V is the seniority quantum number. 
This is relevant only to the nucleir 46 Ti, where for L n — 2 and 4 we have states 
of both seniority 2 and seniority 4. 

The expression for the Ml transition is 



B(Ml)t= (3/4tt)/4 (g P -g n ) x 



]T D° (L, LV) D 1 (L, LV) y/(L(L + 1)) 



LV 



The quantity (g p — g n ) is set equal to 1.89. The coefficients D 1 (L p ,L n ) are 
contained in Ref. 17. The results for the energies and Ml rates are 



E(MeV) 
5.81 
4.00 
4S Ti 3.83 



44 Ti 

46 Tj 



B(M1) f (fj, 2 N ) 
2.70 
1.70 
0.69 



In 44 Ti we have equal numbers of neutrons and protons. All 1 + states in 
the f 7 /2 model have isospin T=l. The J=0 + ground state, of course, has T=0. 
This nucleus is therefore a special case. 

However, for 46 Ti and 48 Ti, the if states have the same isospin as the ground 
state. It is therefore not unreasonable to compare these nuclei with 156 Gd and 
its neighbors. Considering the wide spread in mass number between the nuclei, 
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the behavior is remarkably similar. Perhaps then the presence of 1 + state with 
single particle strength in the 3-4.5 MeV range is a widespread phenomenon. 

Whether the physical interpretation of the 1 + states in the two different 
regions is the same is another question. Independent of what the answer is, 
it will clearly be of interest to try to verify the presence over a wide range of 
the periodic table of 1+ states whose energies and strengths vary systematically 
with mass number. 

The expression for the Ml rate by Dieperink 7 is 

B (Ml) t= ^ (4A^) / ((JV W + N v ) - g v f f N ) 

where N„ and N u are the number of proton and neutron boson, respectively, 
and and g v are the boson g factors. In the single j shell model, the g factors 

are the same for the boson pairs ify 2 ) as f° r ^ ne Sln g ie particle states. 

If we naively apply this formula to the fj/2 shell, with (g n — g u ) set equal 
to (<7p — 9n) = 1-89 then we obtain for 46 Ti. 

B(M1) t= 2.27/4 

and for 48 Ti. 

B(M1) t= 1.71/4 

These are larger than what we calculate in the model. 

In the single j shell model, one cannot ascribe the collective mode as a pure 
orbital mode. This is clear from the fact that we can replace gil + gis by gjj, 
where for j-1+1/2 

9j = (Vj)9i + 9s/ '(2.7) 

while for j-1-1/2 

fc- = [(/ + l)/0' + l)]fli-fl-/[2(j + l)] 

If we use the free values of gi and g s , then the orbital contribution to (g p — g n ) 
is 0.86, and the spin contribution is 1.34. 

The common choices for the renormalized values are gi = 1.1 for a proton 
and -0.1 for a neutron, and g s = 0.7g s (free) for the isovector term. With these 
values, the orbital contribution is 1.03 and the spin part is 0.94. We see that 
the orbit and spin contributions are nearly equal. 

There is one common bond with the IBA formula, 6 ' 7 of course. In both 
expressions one has a factor (g n — g^) 2 . 
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Limited symmetry found by comparing calculated 
magnetic dipole spin and orbital strengths in 4He 
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Department of Physics and Astronomy, Rutgers University Piscataway, New 

Jersey 08855 

(received 17 June 1994) 
Abstract 

Allowing for 2frw admixtures in 4 He we find that the summed magnetic 
dipole isovector orbital and spin strengths are equal. This indicates a 
symmetry which is associated with interchanging the labels of the spin 
with those of the orbit. Where Higher admixtures are included, the orbital 
sum becomes larger than the spin sum, but the sums over the low energy 
region are still nearly the same. 

In Table I we give the total summed strength B(Ml) spin and B(Ml) or b it to 
all (nonspurious) J = 1+, T— 1 states corresponding to the operators ~it z 
and lt z , respectively [as mentioned before we drop the isovector factor 5.586-(- 
3.826)=9.412]. We do this for progressively increasing model spaces: up to 2fko, 
up to Atujj, and up to Qhuj. 

We perform the calculations with the spin-orbit and tensor interactions off 
and on. 

Examining Table I we find one priori unexpected result. When we restrict 
the ground state correlations to 2huj, we find that the summed spin strengths 
are virtually equal to the summed orbital strengths. This is true for all four 
cases of (x,y), i.e., whether or not there is a tensor interaction present. 



Table I. Summed spin and orbital magnetic dipole moment strengths in He 
in units of lCP 3 /z^. 



Inter action 


up to 6huj 


up to Aftw 


up to 2haj 


x y 


Spin Orbit 


Spin Orbit 


Spin Orbit 




1 

1 

1 1 


0.8546 0.8546 
0.8569 0.8571 
3.8245 3.8239 
3.3944 3.3955 


1.3357 5.1635 
1.3417 5.1851 
5.2346 10.937 
4.8288 10.554 


1.5897 7.1474 
1.6211 7.2296 
6.0653 14.607 
5.6052 14.272 



TABLE II. For the cases x=0, y=0 (central interaction, LS limit), we give 
the energies and B(Ml)'s of "spin excited" and "orbit excited" states, with up 
to 2ftu admixtures. 
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Energy (MeV) 


B(Ml) (in 


units of 10 ^ f-^'j\[) 






SptTl 


Orbit 


Nonspurious 


3.67 










44.0 










45.3 


0.855 







48.8 










49.2 










53.9 










56.5 





0.835 


Spurious 


436.7 










436.7 










436.7 










439.3 





13.07 



We consider the case x=y=0. We are in the LS limit. Since the 0s 4 closed 
shell has L=0, S=0, only 2hw excitations with the same quantum numbers will 
admix into the ground state. Let us consider two particles excited from the 
0s shell to the Op shell. We can label the 2p-2h states by[L v L u ] L ~ a [S n S v ] S ~° . 
There are several cases to be considered: 

(1) Two protons are excited. The configurations are (pj) * * (sj) " 
Since L v =0 and S v =0 and L and S are zero, we must have L^={) and SV=0. 
So all in all we get the state \a) = (p 2 ) Ln ~°' " _ ° ^ S ^ L "~ ' "~° _ 

(2) Two neutrons are excited. By analogy, the configuration is | b) = ^ S 2^ L *- > S *- (p 2 ^ L " 

(3) A neutron and a proton are excited from the s shell to the p shell. 
The configuration is {(sp) L7,S7r (sp) LvS "] L=a ' S=a . There are two possibilities: 
\c) = [L„ = \,L V = 1} L=0 [S V = 0,S V = 0} s=0 and \d) = [L n = 1,L„ = 
1] L=0 [S* = 1,S V = 1] S=0 . _^ _^ 

We can form an isovector orbital excitation by applying the operation L v —L v 
to the J=0 + ground state; likewise we can form an isovector spin excitation by 
applying the operator — S v to the J=0 + ground state. When acting on the 
configurations a) or |b), the orbital operator L n — L v gives zero; likewise the 
spin operator — S v . That is, (L n — L V )\L^ = 0, L v = 0) = 0. 

Let us skip to the state |d). Note that the orbital and spin quantum numbers 
are the same L v = S^ = 1 and L v = S v + 1. This is enough to prove that, if this 
were the only state present, we would have the result B(Ml) spin — B(Ml) or bi t t. 

In more detail, (it - lt)\d) = N[L n = l,L v = 1] L=1 [^ = l,S v = l] s=0 
and (St - St)\d) = NIL^ = 1,L v = l] i=0 [^ = 1,S V = l] s=1 . 

There is no reason why these states should be at the same energy and indeed 
they are not, but the equality of the spin and orbital strengths, provided the 
state |c) were not present, is obvious. However, the presence of the state |c) 
apparently presents a problem. The isovector spin operator S n — S v will anni- 
hilate this state, whereas the isovector orbital operator (L n — L v ) creates the 
state [Lk = 1,L u = l] i=1 [5 7r = 0,S U = 0} s=0 . There should therefore be more 
orbital strength than spin strength. What saves the day is that this transition 
is spurious. In the OXBASH program [10] the spurious states are put very 
high in energy by adding a large constant to the single-particle energies for the 
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center of mass motion. We added 100 MeV for each nucleon, thus putting the 
spurious states in the vincinity of 400 MeV excitation energy. In Table III we 
show the 2hw x=0, y=0 calculation in which all the 1+, T=l states are shown, 
both nonspurious and spurious, with the values of B(Ml) spin and B(Ml) or bit- 
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Abstract 

We address some properties of the quadrupole-quadrupole (Q.Q) in- 
teraction in nuclear studies. We here consider how to restore SU(3) sym- 
metry even though we use only coordinate and not momentum terms. 
Using the Hamilton H = Y\ (p 2 /2m + m/2to 2 r 2 ) X J2 l<3 Q W • Q (J) - 
(x/2)Ei<3(*) ' Q(i) with Qm = r 2 Y 2 ^, we find that only 2/3 of the 
single-particle splitting e d — £is comes from the diagonal term of Q.Q; 
the remaining 1/3 comes from the interaction of the valence nucleon with 
the core. The same is true in the Of-lp shell. [S0556-2813(97)01203-X] 

Here we wish to obtain Elliott's SU(3) results [2] in a shell model calculation in 
which only the coordinate Q.Q interaction is used. We do not wish to used the 
momentum-dependent terms. The latter were introduced by Elliott so that, in 
combination with the coordinate terms, there would be no AN=2 admixtures, 
i.e., no admixture from configurations involving 2huj excitations. However, in 
many cases we want to see the effects of such admixtures in our shell model 
studies. One classic problem in which AN=2 admixtures are important is the 
isoscalar E2 effective charge which gets enhanced by a factor of 2 when such 
admixtures are allowed. There are many other problems of interest along these 
lines, some of which we have considered [1]. 
The Hamiltonian we consider is therefore 

H = E(£ + x E g« ■ QU) \ E Q(o ■ g« 

i i<j i 
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where 

Q{i) k . Q(j) k = (-l)V(2fc + l)r(i) k r(j) k [Y(i) k Y(j) k }° 

with k=2. Like Elliott, we have not only the two-body Q.Q term, but also the 
i=j single-particle term. 

It is convenient to introduce the quantity^; = 56 4 x/32-7r where b is the oscil- 
lator length parameter, such that b 2 = h/muj = 41.46/fiw. 

To evaluate the single-particle term we use the addition theorem 

Pk (cos 12 ) = ]T Y k ,» (1) (2) 

and thus 

V5[Y>(i)Y>(i)] = ±-P 2 (l) = ^- 
The single-particle potential is then 

The expectation values of U(r)/x for the single- particle states, Os, Op, Od, 
Is, Of, and lp are, respectively, -15, -35, -63, -75, -99, and -119. What single- 
particle splittingeod — eis is needed to get Elliott's SU(3) results? The best way 
to answer this is to give the formula for the SU(3) energy in the ls-Od shell (in 
which the momentum terms are included): 

E (A/x) = x [-4 (A 2 + f + A/i + 3 (A + M )) + 3L (L + 1)] 

For a rotational band, the L=2-L=0 splitting. 

The splitting due to the diagonal Q.Q interaction is [—63 — (—75)] \ = 12x- 
This is only 2/3 of the desired result. Where does the remaining 1/3 (6%) come 
from? 

The answer is that the missing part comes from the interaction of the valence 
particle with core. That is to say, in order to get Elliott's SU(3) results we must 
not only include the diagonal term, but also the particle-core interaction. 

The expression for the particle-core interaction is 

^ = -xE <*m ( 2 ) \Q ■ Q\ (!) ^ ( 2 ) " < ( 2 ) *°m c (1)) (1) 

c,m c 

where j and m represent all the quantum numbers of the valence nucleon (includ- 
ing isospin labels which have been suppressed) and c and m c are the quantum 
numbers of a particle in the core. 

With the above Q.Q interaction only the exchange term survives. The ex- 
pression becomes 

^ = xE(« c « t ^) (2) 
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We obtain 5eod — S^is = This is missing 1/3 of the splitting required to 
get Elliott's SU(3) results [1]. 

The above results for the ls-Od shell are more general. In the lp-Of shell 
the single-particle splitting required to get the SU(3) result is eo/ — ei p = 
3(3x4=lx2)x = 30x- Once again we only get 2/3 of this 20x from the di- 
agonal Q.Q term. The remaining 10% comes from the interaction of the valence 
nucleon with the core (actually only the Op shell in the core will contribute). 
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Abstract 

For a Q.Q interaction the energy-weighted sum rule for isovector or- 
bital magnetic dipole transitions is proportional to the difference ^B(E2, 
isoscalar)-^B(E2, isovector), not just to YJB(E2, physical). This fact is 
important in ensuring that one gets the correct limit as one goes to nuclei, 
some of which are far from stability, for which one shell (neuron or proton 
is closed. [S0556-2813(97)03208-l] 

PACS number(s): 21.60.Fw, 21.10.Re, 21.60.Ev, 21.60.Cs 

Using the interaction —xQ'Q, Zamick and Zheng [1,2] obtained a sum rule which 
relates the scissors mode excitation rate (i.e., the isovector orbital magnetic 
dipole excitation rate) to the electric quadrupole excitation rate. The isovector 
orbital magnetic dipole operator is (L„ — L v )/2 [the isoscalar one is half the 
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total orbital angular momentum L /2 = (L v 
reads 



— L v )/2\. In more detail, sum rule 



]T (E n - E ) B (Ml) 



■f 16tt 



^ { [B (£2, Oi 2 4 ) /s - B (£72, Oi -+ 2;) 7y ]} (1) 



k 



where B{Ml) n ^ is the value for the isovector orbital Ml operator (g^ = 0.5, gi u = —0.5, g S7I = 0, g sv = 0) 

and the operator for the E2 transitions is£ protons e p r 2 F 2 + £„ eMtrons e n r 2 Y 2 

with e p = 1, e n = 1 for the isoscalar transition(IV). The above result also holds 

if we add a pairing interaction between like particles, i.e., between two neutrons 

and two protons. Our main objectives in this work are to clarify the role of the 

isovector B(E2) in the above formula and to compare the fermion and boson 

model approaches to scissors mode excitations. 

The above work was motivated by the realization from many sources that 
there should be a relation between the scissors mode excitation rate and nuclear 
collectivity. Indeed, the initial picture by Loludice and Palumbo [3] was of an 
excitation in a deformed nucleus in which the symmetry axis of the neutrons vi- 
brated against that of the protons. In 1990-1991 contributions by the Darmstadt 
group [4,5], it was noted that the Sm isotopes, which undergo large changes in 
deformation as a function of mass number, the B(Ml) scissors , was proportional 
to B(Es, 0i — > 2i). The B(E2) in turn is proportional to the square of nuclear 
deformation , 5 2 . 

The above energy-weighted sum rule of Zamick and Zheng [2] was attempt 
to obtain such a relationship microscopically using fermions rather than inter- 
acting bosons. To a large extent they succeeded, but there were some differences 
relative to [4,5]. Rather than being proportional to B(E2, 0i — )■ 2i), the propor- 
tionality factor was the difference in the summed isoscalar and summer isovector 
B(E2)'s. Now one generally expects the isoscalar B(E2), especially to the first 
2+ state, to be the most collective and much larger than the isovector B(E2). If 
the latter is negligible, then indeed one basically has the same relation between 
scissors mode excitations and nuclear collectivity, as empirically observed in the 
Sm isotopes. 

However, the derivation of the above energy-weighted sum rule is quite gen- 
eral and should therefore hold (in the mathematical sense) in all regions, not 
just where the deformation is strong. To best illustrate the need for the isovec- 
tor B(E2), consider a nucleus with a closed shell of neutrons or protons. In such 
a nucleus, and neglecting ground state correlations, the scissors mode excitation 
rate will vanish as one needs both open shell neutrons and protons to get a finite 
scissors mode excitation rate. However, if we have, say, an open shell of protons 
to get a finite scissors mode excitation rate. However, if we have, say an open 
shell of protons and a closed shell of neutrons, the B(E2, 0i — > 2i) can be quite 
substantial. Many vibrational nuclei are of such a type, and they have large 
B(E2)'s from the ground state, e.g., 20 Weisskopf units (W.u.). 

However, in the above circumstances (closed neutron shell), the neutrons 
will not contribute to the B(E2) even if we give them an effective charge. But if 
only the protons contribute, it is clear that B(E2, isovector)=B(E2, isoscalar). 
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As an example, let us consider the even-even Be isotopes 6 Be, 8 Be, 10 Be, and 12 Be. 
In doing so, we go far away from the valley of stability, but this is in line with 
modern interests in radioactive beams. 

BRIEF REPORTS 

TABLE I. The values of B(Ml) orbital , B(E2) isoscalar , B{E2) lsovector for Be 
isotopes. 



Nucleus 




B(fm) a 


B(Ml) or bital 


B ( 1) isoscalar 


B(J\d\ ) isovector 








1.553 





15.63" 


15.63" 


s Be 






1.597 


0.637 


73.54 c 


6.47 


w Be 


T=1^T= 


1 


1.635 


0.0895 


69.7 


31.27 




T=l^T= 


2 




0.149 





3.20 


12 Be 






1.669 





20.85 


20.85 



a b 2 =41.46/(ftw), hw = 45/A 2 / 3 - 25/A 1 / 3 . 

"The analytic expression in 6 Be is B{E2) = (50/47r)6 4 e 2 . 

c The analytic expression in 8 Be is B(E2) = (35/47r)fo 4 (ep + en) 2 . 



Fayache, Sharma, and Zamick [6] have previously considered 8 Be and 10 Be. 
The point was made that these two nuclei had about the same calculated 
B(E2,0\ — > 2i), but the isovector orbital B(Ml)'s were significantly smaller 
in 10 Be than in 8 Be. This went against the systematic that B(Ml)orbital 
is proportional merely to B(E2). In detail, the calculated B(Ml,0i — > 1) 
was {2/tt)h 2 n for 8 Be and in 10 Be was (9/32n)n%(T = 1 -> T = 1) and 
(15/32tt)/x 2 v (T = 1 -> T = 2). Thus the ratio of isovector orbital B(Ml)'s 
is 10 Be/ s Be = 3/8. 
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Abstract 

The g 9 / 2 shell of identical particle is the first one for which one can have 
seniority-mixing effects. We consider three interactions: a delta interac- 
tion that conserves seniority, a quadrupole-quadrupole (Q.Q) interaction 
that does not, and a third one consisting of two-body matrix elements 
taken from experiment ( 98 Cd) that also leads to some seniority mixing. 
We deal with proton holes relative to a Z=50, N=50 core. One surprising 
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result is that, for a four-particle system with total angular momentum 
1=4, there is one state with seniority v=A that is an eigenstate of any 
two-body interaction-seniority conserving or not, The other two states 
are mixtures of v=2 and v=A for the seniority-mixing interactions. The 
same thing holds true for 1=6. Another point of interest is that, in the 
single-j shell approximation, the splitting AE — E(I max ) — E(I m i n ) are 
the same for three and five particles with a seniority conserving interaction 
(a well-known result), but are equal and opposite for a Q.Q interactions. 
The Z=40, N=40 core plus g 9 / 2 neutrons (Zr isotopes) is also considered, 
although it is recognized that the core is deformed. 

III. SPECIAL BEHAVIORS FOR 1=4+ AND 6+ STATES OF THE 
g9/2 CONFIGURATION 

For a system of four identical nucleons in the g 9 / 2 shell, the possible seniorities 
are f =0,2, and 4, with v=Q occurring only for a state of total angular momentum 
1=0. There is also a v=\ state with 1=0. 

For !=4 and 6, we can have three states, one with seniority v=2 and two 
with seniority v=4. For the two v=A states we have at hand, we can construct 
different sets of v=A states by taking linear combinations of the original ones. 
If the original ones are (4)i and (4)2, we can form 

(4) A = a(4)i+6(4) 2) (2) 

(4) B = -6(4)i + a(4) 2 

with a 2 + b 2 = 1. The set (4)A, (4)B is as valid as the original set. 

However, we here note that if we perform a matrix diagonalization with any 
two-body interaction-seniority conserving or not-one state emerges which does 
not depend on what the interaction is. The other two states are, in general, 
mixtures of v=2 and v=A which do depend on the interaction. The values of 
the coefficients of fractional parentage (cfp's) of this unique state of seniority 4 
as shown in Table I. The states of J Q 7^ 4.5 all have seniority v=A state there 
is no admixture of states with J Q = j = 9/2, be they u—1 or v—2>. Again, no 
matter what two-body interaction is used, this 1=4 state remains a unique state. 

TABLE I: A unique J=4, v=4 SFP for j=9/2 



Jo 


(i 3 Joj|j 4 / = 4,i/ = 4) 


1.5 




0.1222 


2.5 




0.0548 


3.5 




0.6170 


4.5 (v = 


1) 


0.0000 


4.5 = 


3) 


0.0000 


5.5 




-0.4043 


6.5 




-0.6148 


7.5 




-0.1597 


8.5 




0.1853 



Amusingly, this state does not appear in the compilation of seniority-classified 
cfp's of Bayman and Lande [20] or de Shalit and Talmi [5]. We should empha- 
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size that, although different, the Bayman-Lande cfp's are perfectly correct (as 
are the ones of de Shalit and Talmi, whose cfp's are also different from those of 
Bayman and Lande [20]). But then, why do they not obtain the unique state 
that we have shown above? Bayman and Lande use group theoretical techniques 
to obtain the cfp's diagonalizing the following Casimir operator for Sp(2j+1): 

G(Sp 2j+1 ) = ^— J2 (-l) fe (2fc+l) 3 /2[[W]° ( 3 ) 

J odd k=l 

where Ug = J2iLi and U is the Racah unit tensor operator 

(*C,\u!;\V m )=5 jj ,(kjqm\j'm') (4) 

The two seniority v=A states are degenerate with such an interaction and, since 
there is no seniority mixing, we can have arbitrary linear combinations of the 4 + 
states. Only by using an interaction which removes the degeneracy and violates 
seniority, do we learn about the special state in Table I. 

IV. THE ENERGY SPLITTING E(I_max )-E(I_min ) WITH A Q.Q 
INTERACTION 

A well-known result for identical particles in a single j shell is that, if one uses a 
seniority-conserving interaction, then the relative spectra of states of the same 
seniority are independent of the number of particles [5-7]. Thus, for n=3 and 
n=5, the seniority v=i states have the same relative spectrum; for n=2, 4, and 
6, the seniority v=2 states have the same spectrum. These results hold, in 
particular, for the delta interaction used here. 

Now the Q.Q interaction does not conserve seniority and the above results 
do not hold. However, we have noticed an interesting result for n=3 and n=5. 
Consider the splitting E{I max ) — E{I m i n }, v=3, where for (79/2? Imax = 

21/2 and 

Imin = 3/2. For a seniority-conserving interaction, AE(n = 5) = AE(n = 3), 
where for a Q.Q interaction, AE(n = 5) = —AE(n = 3). This will be discussed 
quantitatively later. 

X. THE E(I max ) - E(I mm ) SPLITTING FOR n=3 AND n=5: 97 AG 
VERSUS 95 RH AND 83 ZR VERSUS 85 ZR 

As mentioned in a previous section, the splitting AE — E(I max = 21/2+) — 
E(I m in — 3/2+) is the same for three particles as it is for five particles (or three 
holes and five holes) if one has a seniority- conversing interaction. However, for 
a pure Q.Q interaction, we have AE(n = 5) = —AE(n = 3). 
Using the V( 98 Cd) interaction for 97 Ag and 95 Rh, we find 

AE(n = 3) = 0.77058 MeV 

AE(n = 5) = 0.87818 MeV 



21 



They are both positive, an indication that the seniority-conserving delta inter- 
action is much more important than the seniority-violating Q.Q interaction. 

Talmi had previously concluded, from an analysis of hi 1/2 nuclei with a 
closed shell of neutrons (N=82) , that seniority conservation held to a high degree 
[7,31]. 

Unfortunately, for the 59/2 nuclei that we are here considering ( 92 Tc, 95 Rh, 
97 Ag, as well as the zirconium isotopes 83 Zr, 85 Zr, 87 Zr), although the high spin 
states including I = 21/2+ have been identified, the / = 3/2+ states in 97 Ag, 
95 Rh, and 93 Tc, as well as for Zr isotopes. 

For 83 Zr and 85 Zr, with a fitted interaction (despite misgivings of using a 
single j model space), we find for AE = E(I max ) — E(I min ). 

AE( 93 Zr) = 0.48742 MeV 

AE( 95 Zr) = -0.59355 MeV 

They have opposite signs, which shows that for these fitted interactions the 
Q.Q interaction is much more important for this case-neutrons beyond a Z=40, 
N=40 core. 

But it should be emphasized that the / = 3/2 + state is not part of the fit 
because it has not been identified experimentally. If more levels were known in 
the Zr isotopes, and in particular the low spin level / = 3/2+ (but also 5/2+ 
and l/2 + ), then the picture might change. We strongly urge that experimental 
work be done on all the nuclei considered here in order to locate the missing 
states, especially 1 = 3/2^ and also 5/2^. 
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Abstract 

We find a relationship between coefficients of fractional parentage (cfp) 
obtained on the one hand from the principal-parent method and on the 
other hand from a seniority classification. We apply this to the Redmond 
formula which relates n— m+1 cfp's to n-1— >n cfp's where the principal- 
parent classification is used. We transform this to the seniority scheme. 
Our formula differs from the Redmond formula in as much as we have a 
sum over the possible seniorities for the n— m+1 cfp's, whereas Redmond 
has only one term. We show that there are useful applications of both the 
principal-parent and the seniority classification. 

A recursion formula for cfp's due to Redmond [1] is presented in the books of 
de Shalit and Talmi [2] on p. 528, and Talmi [3] on p. 274. It can be written 
as follows 

(n + 1) [j n M)]J\}j n+1 M] J] [r (a 1 J 1 )jJ\}j n+1 [ ao Jo] J] = 
<5 QlQo ^,7o+n(-l) J »+' / V((2Jo + l)(2J 1 + l))E Q2 , 72 J J ] jI 

x [T" 1 (a 2 J 2 )jJ \}j n a J ] [j"" 1 (a 2 J 2 ) j Ji|}j"ai Ji] (7) 

In the above, square bracket designates the principal parent used to calculate 
the cfp. Actually, the principal parent sometimes looses its significance because 
in some cases more than one principal parent can yield the same cfp. In tables 
of cfp's, the principal parent is usually not listed. The quantities in parentheses 
(a J ) are listed. The cfp with (a J ) is the probability amplitude that a system 
of (n+1) identical particles with quantum numbers (a J a ) and a single nucleon. 

2. Relation between principal-parent cfp's and those in the seniority 
scheme We here note a relationship between the overcomplete set of principal- 
parent coefficients of fractional parentage and those with the seniority classifi- 
cation 

[r M)jj\}r +1 M] j] [r M)jj\}j n+1 M] j] = 
E [j n M)jj\}j h+1 m [r (8) 

In the left-hand side above, the first principal parent is formed by adding the 
(n+l)th nucleon to an n-nucleon to an n-nucleon antisymmetric system with 
good seniority and angular momentum (f J ), then coupling the combined 
sytem to a total angular momentum J, and then antisymmetrizing and nor- 
malizing the total wave function. On the right-hand side, the sum over v is a 
sum over all the possible seniorities of the combined (n+1) system and, for a 
given seniority, over all states with that seniority. 

A proof of the above result will be given in Appendix A. 
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